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We consider cold Rydberg atoms in a one-dimensional optical lattice in the Mott regime with a 
single atom per site at zero temperature. An external laser drive with Rabi frequency Q, and laser 
detuning A, creates Rydberg excitations whose dynamics is governed by an effective spin-chain 
model with (quasi) long-range interactions. This system possesses intrinsically a large degree of 
frustration resulting in a ground-state phase diagram in the (A,f2) plane with a rich topology. As 
a function of A, the Rydberg blockade effect gives rise to a series of crystalline phases commen- 
surate with the optical lattice that form a so-called devil's staircase. The Rabi frequency, $1, on 
the other hand, creates quantum fluctuations that eventually lead to a quantum melting of the 
crystalline states. Upon increasing f2, we find that generically a commensurate-incommensurate 
transition to a floating Rydberg crystal occurs first, that supports gapless phonon excitations. For 
even larger f2, dislocations within the floating Rydberg crystal start to proliferate and a second, 
Kosterlitz-Thouless-Nelson-Halperin- Young dislocation-mediated melting transition finally destroys 
the crystalline arrangement of Rydberg excitations. This latter melting transition is generic for one- 
dimensional Rydberg crystals and persists even in the absence of an optical lattice. The floating 
phase and the concomitant transitions can, in principle, be detected by Bragg scattering of light. 

PACS numbers: 05.30.Rt, 32.80.Ee, 37.10.Jk, 64.70.Rh 



I. INTRODUCTION 



The experimental progress in manipulating ultracold 
atomic gases has established these systems by now as 
versatile quantum simulators of basic many-body Hamil- 
tonians with an unprecedented control over coupling 
parameters pi The interactions between ultracold neu- 
tral atoms are, however, usually limited to on-site- or 
contact-interactions due to the s-wave nature of ultra- 
cold collisions. In order to overcome this limitation, there 
has been a considerable effort recently to create atomic 
gases with (quasi) long-range dipolar interactions)^ Gen- 
erally, long-range interactions can give rise to a variety 
of competing ground states and stabilize new quantum 
phases of matter, that are unfamiliar from systems with 
finite-range interactions only. In the context of ultracold 
gases of dipolar atoms various examples have been stud- 
ied theoretically so far, e.g., supersolids,— crystalline 
phasesr'"— liquid crystal and ncmatic phases<^"— 

Apart from using atoms with a static dipole moment, 
long-range forces can also be generated by laser driving 
atoms to dipolar states or, alternatively, to highly ex- 
cited Rydberg states. Such excited Rydberg atoms then 
experience either a dipolar or a van der Waals force, re- 
spectively, that can become quite strong due to their huge 
polarizability that scales as where the principal quan- 
tum number, n, is typically on the order of n ^ 50 in 
current experiments. The interesting aspect of this re- 
alization is that in either case the additional long-range 
interaction acts only between excited atoms. As a con- 
sequence, it becomes energetically unfavorable to excite 
two atoms at short distance - a phenomenon known as 



the dipole-blockade effect.—"— In the presence of many 
atoms, this blockade mechanism gives rise to interest- 
ing collective behavior—"— especially in low-dimcnsional 
systems 1^"— 

In the present work, we consider a regular one- 
dimensional (ID) lattice of Rydberg atoms and study the 
ground-state phase diagram as a function of laser detun- 
ing, A, and Rabi frequency, VI. Interestingly, this ID sys- 
tem possesses intrinsically, as we will argue below, a large 
degree of frustration, and, as a result, the dipole-blockade 
leads here to various liquid-like and crystalline phases 
with the concomitant melting transitions. This system 
was studied in a number of recent publications^^"— i2£ 
that elucidated several aspects of the problem. The den- 
sity of Rydberg excitations and its dependence on the 
Rabi frequency 57 were studied experimentally and theo- 
retically for a trapped gas without an optical lattice^^i^ 
the dynamics of excitations for Rabi frequencies larger 
than the interaction energy was investigated in Rcf. '23, 
and crystalline phases were considered in Refs. i26 - i28i.30l . 
It is the purpose of this work to complement and to com- 
plete the picture of the phase diagram of the system, 
to identify the universality classes of the phase transi- 
tions, and to put it into context with well-established 
results in statistical physics about melting of crystalline 
phases^ii^^ In particular, we point out that as a function 
of laser detuning. A, the Rydberg excitations first con- 
dense into a gapless incommensurate crystalline state in 
agreement with Ref. [26l but we identify the transition to 
be of Berezinskii-Kosterlitz-Thouless^ type in analogy 
to the Kosterlitz-Thouless-Nelson-Halperin- Young sce- 
nario of dislocation-mediated melting of two-dimensional 
crystals on periodic substrates i^i^ As quantum fluctua- 
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tions are further suppressed, for example, by lowering n 
a second, commensurate-incomniensurate (C-IC) transi- 
tion occurs and a Rydberg crystal forms with a period- 
icity that is commensurate with the optical lattice. At 
the two transitions, the character of Rydberg correlations 
changes abruptly which, as we argue, can in principle be 
detected by Bragg scattering of light. In the following, 
we introduce the considered model in some detail and 
discuss the methods applied throughout this work. 



A. Model 

We consider ultracold bosonic atoms in a ID optical 
lattice in the presence of external laser fields that drive a 
transition from the ground-state of an atom to a highly- 
excited Rydberg state. Experimentally, this transition 
occurs via a two-photon process, the intermediate state 
of which can be adiabatically eliminated resulting in an 
effective two-level description in terms of the ground- and 
the excited Rydberg state<^ Thus, we start from a Bose- 
Hubbard model of atoms in a ID optical lattice together 
with an additional external driving field at effective Rabi 
frequency Q and detuning A. We assign to each atom 
a pseudospin a =],, f that indicates whether it is in the 
ground- or excited Rydberg state, respectively. The laser 
excites Rydberg atoms, and excited atoms can decay by 
spontaneous emission or colhsional processes giving rise 
to a decay rate 1 /r. In general, it is therefore not possible 
to describe the atomic subsystem by a Hamiltonian but 
a description in terms of a density operator must be used 
instead. However, if we limit ourselves to time-scales 
much smaller than the life-time t of the excited Rydberg 
atoms we can approximate the evolution of the density 
operator with the help of an effective Hamiltonian, whose 
ground state phase diagram can be analyzed. 

Employing the rotating wave approximation^ this 
Hamiltonian takes the form 

+ ^n"jT"ji + ^ H -^ii-ei ^jfrin^ (1) 

where hj^ is the bosonic annihilation operator at lattice 
site j and nj^ = b^j^bj^ is the density operator for ground- 
and excited-state atoms at site j. Generally, the hopping 
of atoms, tcr, and their on-site interaction, Uaa' , depend 
on their internal energy state a. In addition, the excited 
atoms experience a repulsive interaction denoted by J|r| 
that decays algebraically with dimcnsionless distance r 



where the exponent is either a = 6 if the interaction is of 
van der Waals type or a = 3 for a dipolar interaction be- 
tween excited Rydberg atoms.— Although the interaction 
of Eq. (121) decays rather fast for the exponents of interest, 
we call it nevertheless (quasi) long-range in this work in 
the sense that it is finite for all distances r. 

We assume that the filling is exactly one atom per site, 
and we consider the Mott regime where the strong on- 
site U ^ to- prohibits the atoms to hop and freezes them 
to the lattice positions. Projecting the "high-energy" 
Hamiltonian ([1]) onto the low-energy Hilbert space with a 
fixed density of one atom per site one obtains an effective 
spin Hamiltonian 

n=Y,J\i^,\{S^ + 1/2) {S^ + 1/2) (3) 

e>j 

- ^ E i^ts^, + h.c.) -Aj2s^ + nY,s^, 

with S*" denoting spin-1/2 operators at lattice site j. The 
Rabi frequency $7 > corresponds to a transversal mag- 
netic field and the detuning A is an additional magnetic 
field component in longitudinal direction. Virtual hop- 
ping processes in the Mott regime give rise to a transverse 
spin interaction Jl ^ /U between nearest neighbors; 
see appendix [Xj Note that J_l > is ferromagnetic as 
we consider bosonic atoms.— In appendix |^ we also 
shortly discuss the opposite limit of a Bosc-Einstcin con- 
densate of ground-state atoms and show that the effective 
Hamiltonian in this case is identified with a spin model 
similar to Eq. as well. 

Eq. is the Hamiltonian that we analyze in this work. 
In the following, we will use the language of spins and 
Rydberg atoms interchangeably. For example, the mag- 
netization per site plus one-half, < to-|- 1/2 < 1, corre- 
sponds to the mean-number of Rydberg excitations per 
lattice site. We will be mainly interested in the ground 
state phase diagram of the Hamiltonian ([3]) in the (A, Vl)- 
plane. Note that this phase diagram is symmetric with 
respect to the line A = A = X]r>o '^hich corresponds 
to vanishing longitudinal magnetic field. In the follow- 
ing we restrict, however, our discussion to the relevant 
regime |A| < Jr, i.e., |A| < A. 

In usual experimental settings the Rydberg interaction 
Jr is by far the largest energy scale, 

C!,|A|,Ji« Jij. (4) 

Previous studies have been limited to the so-called frozen 
Rydberg gas limit where the hopping of Rydberg excita- 
tions Jj^ is omitted arguing that relevant experimental 
time-scales like the Rydberg life-time r are usually much 
smaller than h/J±. Although the bare value of J± is 
indeed small, it becomes strongly enhanced by renormal- 
ization effects generated in second order in fii^ We there- 
fore find it convenient to keep it explicitly in the model 
as it simplifies the theoretical analysis. 
The analysis of Eq. ([3]) in the parameter regime (|1]) is 
actually a formidable task for the following reason. As a 
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starting point of an analysis, we might neglect the fields, 
r2 and A, the transverse interaction Jl and also all lon- 
gitudinal interactions between neighbors with a distance 
larger than a single lattice spacing because of the strong 
decay of J^^^/Jfj = l/\r\"' with increasing distance r. In 
this case, the Hamiltonian reduces to 



(5) 



It turns out that the ground state of Eq. ([5]) is macroscop- 
ically degenerate with zero-point entropy S = log ^-^-^ , 
given by the logarithm of the golden meau;^ as all states 
without two adjacent up-spins, i.e., two adjacent Ryd- 
berg excitations, have the same energy. Almost any ad- 
ditional term will immediately quench this entropy, and 
it is clear that a perturbative approach in the parameters 
A, ri, J_L and J\r\ with r > 1 is, therefore, inadequate to 
derive the properties of the Hamiltonian ([3]). 



B. Methods 

One can obtain an idea of the complexity of Eq. ([3]) 
by considering the classical limit = J± = 0. The 
ground states of the resulting classical model were deter- 
mined by Bak and Bruinsma^ who demonstrated that 
the phase diagram as a function of A possesses infinitely 
many phases, and the magnetization m(A) forms a com- 
plete devil's staircase. For A < the ground state is 
fully polarized ... ■ll-l-l ... with magnetization m = 
per site. For positive A, it becomes energetically favor- 
able to accommodate up-spins, i.e., Rydberg excitations 
in the ground state. The up-spins are however arranged 
in a way so that their distance is maximized thus mini- 
mizing their interaction energy. The resulting mean dis- 
tance between two Rydberg excitations, known as the 
Rydberg blockade radius, depends on the density of ex- 
citations. By increasing A, more Rydberg excitations 
are accomodated, the blockade radius shrinks and the 
magnetization increases in steps giving rise to a plethora 
of phases that are commensurate with the periodicity of 
the optical lattice potential. The magnetization finally 
vanishes at sufficiently large detuning A for the antifer- 
romagnetic phase ... ititit •■•j where every second atom 
is excited. Generally, the commensurate phases break a 
certain Zp symmetry of the optical lattice where p is the 
periodicity of the spin configuration. For example, the 
antiferromagnetic state repeats itself after every second 
site so that it breaks an Ising, i.e., a Z2 symmetry. The 
devilish aspect of the phase diagram is due to its fractal 
appearance as phases of all possible commensurabilities 
with the optical lattice are realized on the classical level. 

The presence of a finite and J± induce 
quantum fiuctuations, which eventually melt the 
commensurate phases, and one enters the realm 
of commensurate-incommensurate transitions, floating 
phases and dislocation-mediated melting intensively 



studied at the dawn of the 1980's,^ii^ We are going to ap- 
ply the results and techniques developed at that time to 
the present context. In particular, throughout this work 
we apply intensively the strategy of Villain and Bak^ as 
used in their analysis of the axial next-nearest-neighbor 
Ising (ANNNI) model. First, we consider the classical 
Hamiltonian close to one of its phase transitions and an- 
alyze the effect of a finite hopping Jj_ perturbatively still 
keeping f2 = 0. This allows to derive an effective low- 
energy Hamiltonian for excitations of the Rydberg crys- 
tal - domain walls - that represents a proper fixed-point 
theory in the renormalization group sense. We are then 
able to apply a standard stability analysis of this effec- 
tive theory with respect to quantum fluctuations induced 
by a finite Rabi frequency f2, that creates topological de- 
fects, i.e., dislocations. This strategy allows us to derive 
the phase diagram in the limit f2 <C J_l. We then com- 
bine this knowledge with known results in order to derive 
the topology of the full phase diagram of the Hamilto- 
nian Eq. in particular, in the experimentally relevant 
regime V,:^ J±. 

The rest of the article is organized as follows. In section 
|TI]we consider the various gapped commensurate phases 
where the Rydberg crystal is commensurate with the op- 
tical lattice and construct an exact critical theory for the 
transitions between those crystalline states. In section HIIl 
we consider the Hamiltonian ([3]) in the continuum limit 
and discuss the melting of the incommensurate floating 
phase. We conclude with a summary and discussion in 
section IIVI In appendix |X] we shortly discuss the su- 
perfluid regime, which has been addressed in a recent 
experiment on Rydberg atoms in a ID latticci^ In ap- 
pendix [B] wc establish the effective particle formulation 
used in the text to describe arbitrary transitions between 
classical commensurate phases. 



II. COMMENSURATE RYDBERG-CRYSTAL 
PHASES 

In this section, we derive the topology of the phase 
diagram of the Hamiltonian in the (A, Q) plane with 
an emphasis on the various Rydberg crystal phases that 
are commensurate with the optical lattice. 

In order to approach the problem and to circumvent 
the complexity of the Hamiltonian ^ at its classical level 
= J_L = 0, where it possesses a fractal phase diagram)^ 
it is useful to consider first a class of auxiliary Hamilto- 
nians of the same form as Eq. ([3]) but the infinite-range 
interaction, J|,.|. of Eq. ^ replaced by an interaction 



J^^^ of finite range only 



jI^i > if l<\r\<n, 
Jf;^ = if |r| > n. 



(6) 



This finite range interaction should have the property 

(n) (n) 

that it obeys the convexity condition, J^r-i\ + "^|V-i-i| ^ 
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2J| J| for each n and that it coincides with J|,|,| in the 
hmit ri — >■ oo. We start our analysis by considering only 
interactions between close neighbors with small ri, and 
afterwards we discuss the modifications obtained by in- 
cluding interaction of higher order by increasing n. This 
allows us to develop the topology of the phase diagram 
step by step. The full Hamiltonian is obtained in the 
limit n — > (X), but it turns out that for a finite hopping 
Jj_ the topology of the phase diagram in the (A, il) plane 
becomes insensitive to interaction components J|r| at a 
sufficiently large distance r so that J|,,| <^ J±. 

The case n = I when the Rydberg interaction is lim- 
ited to nearest-neighbors only, J|^^"'^^'' = jRS\r\,i, is spe- 
cial and well-known, so we summarize here only the re- 
sult. In this case, the Hamiltonian reduces to a standard 
strongly anisotropic spin-chain in a strong longitudinal 
and a weak transversal magnetic field. In the classical 
limit O = Jl = 0, there are just two ground states, the 
antiferromagnetic state ... titit ... for A > and the 
fully polarized state ... ... for A < 0. At the tran- 

sition A = 0, the classical Hamiltonian reduces to Eq. (O 
with a macroscopically degenerate ground state manifold. 
The corresponding zero-point entropy, however, is imme- 
diately quenched by an infinitesimal hopping Jl and the 
single transition is replaced by an extended Luttinger 
liquid phase bounded by two Lifshitz transitions^^ The 
Luttinger liquid phase, however, is itself unstable with 
respect to a finite transverse field f2, and it gives way to 
a single line of Ising transitions in the (51, A) plane4^ 

This picture, however, changes qualitatively upon in- 
cluding Rydberg interactions beyond nearest neighbors. 
In the next section Hi A[ we consider the spin-chain with 
nearest and next-nearest neighbor interactions before dis- 
cussing the long-range interaction in section IIIBI 



A. Anisotropic spin-chain with nearest and 
next-nearest neighbor interaction 

The auxiliary Hamiltonian with the infinite range in- 
teraction replaced by J|^| of Eq. ©, i.e., with nearest, 
Jl, and next- nearest neighbor interaction, J2, reads 

n=Y.T.Jr{S] + 1/2) + 1/2) (7) 

j r = 1.2 

3 3 3 

It will be important that the interactions obey the con- 
dition of convexity Ji > 2J2 > 0; see Eq. ([2]). In this 
section, we discuss in detail the phase diagram of Eq. ([7]) 
in the limit $7, |A|, Jj^ ^ Ji. We start with a discus- 
sion of the classical ground states at $7 = Jj_ = and 
identify two phase transitions. We then analyze how the 
phase transitions are influenced by a flnite J± and Q and 
present the resulting phase diagram. 



The discussion of this section is closely related to work 
by Fendley, Sengupta and SachdevJ^ In fact, the Hamil- 
tonian d?]) in the limit J_l — )■ and Ji — )• 00 reduces to 
the one studied in Ref. However, the presence of a 
flnite Jl in our case facilitates some of the analysis and 
allows to present the physics in a transparent manner. 

1. Classical analysis for = J± — 

At ^l = J± ^ 0, there are three classical ground states 
as a function of the longitudinal field A that are sepa- 
rated by two classical phase transitions located at 

Ael = 0, Ae2 = 3J2. (8) 

The classical ground states are the fully polarized state 
.. .. at A < Aci with magnetization m = —1/2, the 

state .. .. with periodicity p = 3 and magneti- 

zation m — —1/6 for intermediate fields Ad < A < Ac2, 
and the antiferromagnet .. ^.f^tit •■ with periodicity 
p = 2 and magnetization m = for A > Ac2. The 
two states with period p break a Zp symmetry of the op- 
tical lattice; the state with m = — 1/6 is a Z3 symmetry- 
broken state and the antiferromagnetic state breaks an 
Ising Z2 symmetry. 

At both transitions, there is a finite residual entropy at 
zero temperature. For example, the lowest energy excita- 
tion of the antiferromagnet is a single pair of two adjacent 
down spins, i.e., a domain wall or a spinon .. titiitit ■•• 
The excitation energy of these domain walls vanishes at 
the second critical field Ac2 so that they proliferate, re- 
sulting in a large degeneracy of the ground state at this 
critical point. The domain walls are strongly interacting 
as they must be separated at least by a single down spin. 
The corresponding entropy per site can be evaluated by 
the transfer matrix method that gives in the thermody- 
namic limit L — > 00 

5 = iloga; = 0.28119... at A = Ac2, (9) 

where x is the real root of the cubic equation — 1 -I- .t — 
2x^ + a;'^ = 0. Similarly, the lowest energy excitation of 
the fully polarized state is a single up spin .. mfiii •■ 
whose excitation energy vanishes at Ad =0. At the 
critical field Ad all states with up spins separated at least 
by two down spins are degenerate leading to a residual 
entropy 

= i log a; = 0.38224... at A = Ad = 0, (10) 

where x is the real root of the equation —l — 2x—x^+x^ = 
0. 

In the presence of a finite J± the excitations acquire 
kinetic energy and each single classical transition is re- 
placed by two C-IC transitions enclosing an extended 
floating Luttinger liquid phase; see Fig. [T] We flrst dis- 
cuss these transitions close to the second critical field Ac2 
and afterwards the ones close to Ad. 
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K = 1\ /K = i A' = 

.. iitiiTiiT 




A,i = 



FIG. 1: Phase diagram in the (A, Jx) plane of the auxil- 
iary Hamiltonian ([7| at Q = 0. A finite J± gives rise to 
Luttinger liquid phases (shaded areas) emanating from the 
classical transition points at Ad — and Ac2 = 3J2. The 
spin configurations indicate the classical ground states in the 
three phases with magnetizations m = —1/2, m = —1/6 and 
m = per site, respectively. The values of the Luttinger pa- 
rameter K at the four Lifshitz transitions given in the figure 
is derived in the text. 



2. Transition at Ac2 

We first focus on the transition at Ac2 and study the 
influence of a finite hopping J± that leads to an extended 
Luttinger liquid phase. Afterwards, the stability of these 
phases is analyzed with respect to a finite driving fi. 

a. effective Hamiltonian for small Jj^ > and 17 = 0. 
At the classical transition located at Ac2 the ground state 
changes its period from p = 2 to p = 3. Generally, the 
unit cell of a commensurate state with periodicity p in- 
volves p spins that we denote as Up, for example Up —\.'\ 
for p = 2 and Upi =i^t for = 3. Near the classical tran- 
sition between two adjacent commensurate phases we can 
limit ourselves to the degenerate states with the succes- 
sive occurrence of unit cells of the p = 2 and = 3 states 
only, e.g. ..U2U2U3U2U3U3... The border between two dif- 
ferent unit cells corresponds to a domain wall. In order to 
capture the meandering of domain walls we follow Villain 
and Bak^i^ and introduce fictitious particles living on a 
fictitious lattice to describe the commensurability transi- 
tion. We associate the unit cell U2 with an empty site and 
the unit cell M3 with an occupied site so that, e.g., the 
above sequence of unit cells is identified with ...001011... 
The number of particles is denoted as N and the number 
of empty sites by so that L = 2Ne + 3N where L is the 
length of the chain. The number of fictitious particles is 
bounded by < A'' < i/3 and the length of the fictitious 
lattice is L f = Ne + N = {L — N)/2. The magnetization 
per site is a function of A^ and given by m = —N/ {2L) 
ranging between — 1/6 < jti < 0. 

At the classical transition, all states with < A^ < -^/3 
are degenerate. We can evaluate the resulting entropy 
S = log fis where fi^ is the total number of available 
states. 



exactly reproduces Eq. ([9]) in the limit L 00. Note 
that we are ignoring here boundary effects that depend 
on the employed boundary conditions and arc negligible 
in the thermodynamic limit. 

A finite J± allows the particles to propagate on the 
fictitious lattice. The effective Hamiltonian describing 
this propagation of domain walls is, in lowest order in 
J±, just the tight-binding model for hard-core bosons 



i=l L 



A-A,2 t 



(12) 



with the hopping amplitude J_l/2. All states must obey 
the hard-core constraint ala.^ < 1. The effective chem- 
ical potential is obtained by realizing that each particle 
carries the magnetization — i and, in addition, that at 
the critical field A = Ac2 for J± = the energy is inde- 
pendent of the number of particles A^. 

The ground state energy in the presence of A^ particles 
is easily evaluated. 



,.A-A^ ^ f2Tr 

E = N Jj^ > cos — n 

n=-N/2 ^ ^ ■ 



(13) 



Using that A^ = —2Lm one obtains the magnetic energy 
per site, e = E/ L, as a function of m that reads in the 
limit L — > CX3 



/ N / . AN Ja- , ^ . f 47rm 

e m = - m A - Ae2 + —(1 + 2m sm — 

27r V 1 + 2m 



(14) 



Minimizing e{m) with respect to the magnetization 
— 1/6 < m < 0, one obtains the magnetic equation of 
state de/dm ~ 0. From the second derivative, one gets 
the differential susceptibility = ff^e/dm?, 



X = 



(1 + 2mf 



SttJ^ sin 



\ l+2m ) 



(15) 



The classical jump of the magnetization from zero mag- 
netization to m — —1/6, implying an infinite susceptibil- 
ity at the classical transition, is smoothened out by the 
finite hopping Jj_. Instead, the classical transition is re- 



placed by two Lifshitz transitions at the fields A 
with 



A 



L/3 



^ ATI (k^ 

Af=0 V 2 



(11) 



A+(f7 = 0) = A,2 + 2Jx, 
A-2(l] = 0) = A,2-3Ji. 



(16) 
(17) 



where each term in the sum is just the number of ar- 
rangements of A^ particles on Ly = {L — N)/2 sites. This 



Close to these critical fields the magnetization 
changes with the square-root cusps typical for C- 
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IC transitionsi^i^ 



47r V Ji_ 



J± 



(18) 



6 97r V 3 J± 



(19) 



For fields A~2 < A < A^, the gas of domain walls is 
compressible and forms a Luttinger liquid, whose prop- 
erties we consider in the following. 

b. Luttinger liquid theory for J± > and its stability 
with respect to Q perturbations. The low-energy prop- 
erties of the effective Hamiltonian ([T^ in the magnetic 



field range A^2 
liquid theory 



< A < A^ is captured by the Luttinger 



Hll = 1^ y dx (^Kidj f + -|(a,0)') , 



(20) 



where the bosonic field represents fluctuations in the 
magnetization, 6m ^ —^dx(f>, and [^dx4>{x),9{x')] = 
~iS{x — x'); the two parameters i; and K are to be de- 
termined. The velocity v is given by the energy spacing 
for a finite size system Al = ^v. From the Hamiltonian 
([T2|) . we get Al = ^TJ' sm{Nn/Lf) and after expressing 
N and L f in terms of the magnetization we obtain 



2Jj^ . -Airm 
■ sm ■ 



1 + 2m 



1 -t- 2m 



(21) 



From the magnetic susceptibility (|T5|) and the relation^ 
X = K/(t:v) we obtain the important result for the Lut- 
tinger parameter 



(22) 



As the upper boundary A^2 of the Luttinger liquid re- 
gion is approached, the magnetization vanishes and the 
Luttinger parameter assumes the value K = 1/4. Close 
to the lower boundary A~2 the magnetization becomes 
m = —1/6 and thus K = 1/9; see also Fig. [1] The val- 
ues of the Luttinger parameter close to AJ2 consistent 
with the general expectations for Lifshitz transitions into 
commensurate states of periodicity p that are stabilized 
by Umklapp scattering processes of p spinous yielding 
K = l/p'^ 

Physically, the Luttinger liquid phase corresponds to 
an incommensurate, floating crystalline state of Rydbcrg 
excitations. It has no true long range order but alge- 
braically decaying crystalline correlations^ 



{S'{x)S'{Q)) -m" 



cos [7r(l -I- 2m)x] 



\2K 



(23) 



Now we are in a position to analyze the effect of the 
transversal magnetic field f2. Within the Luttinger liq- 
uid framework, the transversal field gives rise to the 
perturbation^ 



n 



pert 



27ra 



dx cos(6'). 



(24) 



where a is the optical lattice spacing. Note that the inter- 
action J_L is ferromagnetic so that the low-energy fluctu- 
ations of the operator are non-staggered. The origin 
of this perturbation becomes apparent by realizing that 
the conjugate fields 4){x) and 9{x) can be viewed as the 
polar- and azimuthal angles of classical vectors represent- 
ing the spins. The transverse field term ((24)) thus wants 
to lock the azimuthal angles at = or tt, depending 
on the sign of f2, such that all spins point in x-direction, 
as expected. Using standard scaling arguments;^ one 
finds that the perturbation (|24)) is relevant in the renor- 
malization group (RG) sense if the Lutttinger parameter 
exceeds a critical value K > 1/8. Using Eq. this im- 
plies for the magnetization m < — (2 — •\/2)/4 and finally 
A > A^""" where the critical field is determined by the 
magnetic equation of state 



A^T ^ Ae2 - J± _ 2x/2cos(2y2^) j 



(25) 



Ac2 - Jj_ X 2.5908. 



So the Luttinger liquid generated by Jl is only stable 
very close to the first Lifshitz transition at A^2 i^i^ ^ small 
field range obeying -3 < (A - Ac2)/J_l < -2.5908... 
At the critical field A^"'' the operator ([24| leads to a 
Kosterlitz-Thouless transition into an incommensurate 
gapped state. This gapped phase extend to finite $7, and 
the phase boundary in the (A, il) plane starts at A^"^ 
linearly with $7, which follows from the separatrix of the 
Kosterlitz-Thouless RG fiow^ 

c. Effective Hamiltonian for small finite J± and D,. 
It is instructive to consider the modification of the ef- 
fective model (fT2|) due to the transversal magnetic field 
explicitly with the help of a Schrieffer-Wolf transforma- 
tion. Including corrections up to first order in J± and 
second order in ft, i.e., 0{J±, fP), the low-energy Hamil- 
tonian assumes the form 



'HcS — 



(26) 



where we remind that the size of the fictitious lattice 
Lf = Ne + N = (L — N)/2 depends on the number 
of particles N = X^i'^I'^i- The hopping amplitude t at 
A = Ac2 is given by 



t 



1 



2 V 2 / V Jl - 2 J2 2 J2 



(27) 
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To obtain the chemical potential /i and the interaction 
U, it is useful to consider the correction to the energy 
per site, Se ~ SE/L, of various specific states in order 
0{n^, A — Ac2) in the absence of any hopping t, 



de 



0000.. 



/ 1 



f 



<5£..iiii.. = 



2 I 2 \2J1-3J2 J2J ' 



Se 



6 

A- A,2 



2 J 3 V Ji - 2 J2 3 J; 



010101.. 



10 

nV 1 f 1 



(28) 

)• 

(29) 
(30) 



2 J 5 \ 2Ji — SJ2 J2 Ji — 2J2 



Equation ((28)) gives the correction to the vacuum with 
no particles ..0000.. = ..U2U2U2M2-- where U2 =it- The 
energy correction fe.. 010101.. for the state with particles 
at every second site can be used to compute the chemical 
potential; for ~ N which implies N ~ L/5 wc have 
the relation —fj,N = ((5e.. 010101.. ~ ^..0000..)^- With this 
we obtain for the chemical potential 



Ac2 - A 



5(fc. 



Se 



0000.. ~ "fc.. 010101 



.) 



(31) 



Ac2 - A n'^ 



3 4^3 
2Ji — SJ2 Ji — 2J2 J2 



The interaction U can be obtained by relating the energy 
corrections to the vacuum and to the state with particles 
at every site, <5£..iiii..L = Se..omw..L + N{—ii + U) where 
the number of particles is now maximal, i.e., N = L/3. 
This yields an attractive interaction given by 



U = fi + 3(fe..iiii.. — Se, 



0000, 



.) = - 



Alternatively, we can consider the effective Hamilto- 
nian after a particle-hole transformation where holes are 
created by hi within the fully occupied chain ..111111.., 



(33) 



Again, all states must obey the hard-core constraint 



hlh, < 1. 



The hopping amplitude is the same 



as for the particles = t. However, the effec- 
tive chemical potential, fie, and the interaction, Ug, 
for the holes are now obtained from the relations 
-fieNe = (fe.. 010101.. - '5e..iiii..)5A^e and (Se..oooo..2iVe = 
fc. . 11111. .2A''e-l-A^e(—Me + C^e), respectively. In particular, 
this yields the relation = — 2/x/3 -I- bU/i. 

At this order of perturbation theory, the model ([25)1 
still exhibits two Lifshitz transitions. To obtain the upper 



critical field we can neglect the interaction and simply set 
2t -I- ^ = 0, yielding 



A+(0) = Ae2 +2,/j_ 



(34) 



+ '2 



1 



2Ji — 3J2 Ji — 2J2 J2 



The Lifshitz transition at the lower critical field takes 
place if 2te + /-te = 0, giving 



A;2(^^) = Ae2-3J± 

2 



(35) 
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2Ji — 3o^2 Ji — 2,J2 3J2 



However, the picture changes qualitatively if correc- 
tions of order 0{J±il.) or of order 0{n^) are included as 
processes are generated in the effective description that 
do not conserve the number of particles. For example, 
applying three spin flips two effective particles can be 
created out of three adjacent holes, so that the state 
.. iUUt ■■ = --OOO.. is converted to .. UtUt ■• = •■H--- 
In the dense limit close to A~2 where there are only few 
holes in the ground-state, this process corresponds to the 
additional operator hihi^ihi^2 in the Hamiltonian that 
destroys three holes. This operator is irrelevant and the 
Lifshitz transition near A~2 remains unaffected. In the 
dilute limit, however, this process corresponds to the ad- 
ditional operator alal_^-^ that creates two particles which 
is relevant and converts the Lifshitz transition at Ajtj into 
an Ising transition. This is analogous to the mechanism 
at play close to the multicritical point in the ANNNI 
model where also an Ising transition emerges4^ As a re- 
sult, there is also an additional transition of Kosterlitz- 
Thouless type at Af^, with < Af^ < A+2, be- 
tween a gapped ground state and the Luttingcr liquid 
as discussed in Section fll A 2 bl In order to evaluate the 
correction to A^""" in second order in ft, however, the 
interacting Hamiltonian ([25)1 must be solved4^ 



3. Transition at Ad 

We now turn to the transition at Ad between the fully 
polarized state and the commensurate phase with mag- 
netization m = — i; see Fig.[TJ Again, we first study the 
influence of a finite hopping Jj_ and afterwards we per- 
form a stability analysis with respect to f2 perturbations. 

a. Effective Hamiltonian for small J± > and CI ~ 0. 
We can repeat the reasoning of the previous section 
and derive an effective Hamiltonian also close to the 
transition at Ad = 0. Whereas the ground state for 
Ad < A < Ac2 has periodicity p = 3 with unit cell 
U3 =iit, the ground-state for A < is fully polar- 
ized. The unit cell of the latter is trivial as it consists 
of a single spin only, ui At the classical transi- 

tion A = for Jj^ = 51 = 0, there is again a class 
of degenerate ground states that correspond to the suc- 
cessive occurrence of the two unit cells W3 and ui, e.g.. 
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..U3U3U1U3U1U1.. We introduce effective particles and as- 
sociate Ml with an empty and with an occupied site so 
that the above sequence becomes ..110100... The length 
of the chain is related to the number of particles, N, 
with < < L/3, and the number of empty sites by 
L = Ne + 3N. The particles live on a fictitious lattice of 
size Lf = Ne + N ~ L — 2N and the magnetization per 
site is given by m ^~{Ne + N)/{2L) = N/L -1/2 that 
varies between —1/2 < m < —1/6. 

In the classical limit, there is again a macroscopic de- 
generacy at the critical field A^i with a finite entropy per 
site S = logfis, 



(L-2jV)! 

m{L-2Ny.' 



(36) 



where each term corresponds to the number of arrange- 
ment of N particles on Lf sites. This reproduces the 
entropy of Eq. ([TU]) . 

A finite J± results in hopping of particles described by 
the tight-binding model of hard-core bosons. 



i=i 



J± 
2 



a] a. 



i+i 



- (A - Aci)ala, 



(37) 



where the chemical potential is obtained with the same 
reasoning as for Eq. (|12l) . Repeating the steps of the pre- 
vious analysis in section 111 A 21 we obtain for the mag- 
netic energy per site 



e(m) = — to(A — Aci) 



2J_L 
TT 



'""'(id 



(38) 



Minimization with respect to m gives the magnetic equa- 
tion of state, de/dm = 0, that relates A to the magne- 
tization m in the range —1/2 < m < —1/6. For the 
susceptibility = d'^e/drn?, we obtain 



•^J-i-cosh^) 



(39) 



The hopping J± smoothens out the classical transition 
and gives rise to two Lifshitz transitions located at 

A-i(f7 = 0) = -Jj_, and A+ (O = 0) = 3Jj_, (40) 

close to which the magnetization again varies as a square 
root of A — A~^ and A^^ — A, respectively. 

b. Luttinger liquid theory for Jj^ > and its stabil- 
ity with respect to Q. perturbations. Between the critical 
fields A~^ < A < A^]^, the magnetization is compress- 
ible and its fiuctuations are described by a Luttinger liq- 
uid; see Eq. ([20]) . In the following, we determine the 
velocity v and Luttinger parameter K . The velocity is 
again obtained from the finite size spectrum of Eq. ([57| , 
Al = 2ttv/L, giving 



J_L / TT \ 

- — cos - — 
2m V4m/ 



(41) 



With the help of the relation for the susceptibility^ x = 
K / {■nv) we obtain the important result for the Luttinger 
parameter 



K = {2mf 



(42) 



The Luttinger parameter varies within the compressible 
regime between K = 1/9 close to the upper critical field 
A^j^ where m = —1/6, and K = 1 close to A~^ where the 
magnetization is saturated, to = —1/2; see Fig. [TJ 

A finite transversal field $7 gives rise to the operator 
([24| that becomes relevant for K > 1/8 implying for the 
magnetization the bound m > — 1/ \/32 for the stability 
of the Luttinger liquid. Using the equation of state, we 
obtain the corresponding bound for the field A > A 
with 



KT 



^KT = - ^ (cos(\/27r) + \/27rsin(\/27r)) (43) 
= J_L X 2.89583... 

The Luttinger liquid is only stable within a very narrow 
range 2.8958.. < A/ J_l < 3. At the critical field A^^^ ^he 
operator ([24|) results in a Kosterlitz-Thouless transition 
destroying the Luttinger liquid. 

c. Effective Hamiltonian close to A^j^ for small 
Q.,J± > Q. In the following, we discuss correction to 
the effective Hamiltonian ([37| that are generated by a 
finite transversal field f2. In the presence of f2, the num- 
ber of fictitious particles is not conserved anymore al- 
ready in linear order in 51. In the dilute limit close to 
A~^, the transversal field VL can create particles by flip- 
ping a single spin only, i.e., it can produce M3 unit cells 
within the fully polarized state .. iiii ■■ ~ ..uiuiuiui... 
In first order in J7, an operator of the form a\ is gen- 
erated in the Hamiltonian that is relevant and immedi- 
ately destabilizes the Luttinger liquid in agreement with 
Eq. (j43|) . In the dense limit close to A^^, on the other 
hand, the transversal field can convert in linear order in 

a particle, i.e., a M3 =i4,t unit cell by down-flipping 
the last spin into three holes. This process corresponds 
to an additional term in the Hamiltonian that consists of 
the product of three hole-creation operators h\h\j^^h\j^2- 
This operator is irrelevant in the limit of small densities 
of holes, and the Luttinger liquid survives in a finite sliver 
close to A^^. 

This actually allows to determine the f2-dependence 
of the phase boundary A^j^ in the limit f2 — 0. It can 
be determined by solving the Schrodinger equation for a 
single hole in the background of particles in lowest or- 
der in VL. The excitation energy of the hole vanishes at 
A^]^. Equivalently, we can consider the effective single 
hole Hamiltonian that captures the Lifshitz transition 



(44) 



The three hole-creation operator generated in linear or- 
der in fl and the interaction between holes can be ne- 
glected as they are irrelevant in the RG sense close to 
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KT 




FIG. 2: Sketch of the phase diagram m the (A, SI) plane of 
the auxihary Hamiltonian ((7|) at a smaU but finite hopping 
Jx with two symmetry broken phases (yellow shaded). A 
broken Z2 phase is bounded by an Ising transition (/) and 
a broken phase is enclosed by two Lifshitz (L), i.e., C-IC 
transitions except at a single point where the transition is 
in the Z-i Potts universality class (P3). The Luttinger liquid 
phases (blue shaded) are confined to tiny regions around the 
Z3 phase and are separated from the incommensurate gapped 
phase by a Kosterlitz-Thouless transition (KT). 



the Lifshitz transition. The hopping and the chemical 
potential in lowest order in f2 can be evaluated similar as 
in the previous section, 



j_L fnyfi 1 



2 V2/ V^ -^i-^ 
A fQ.V f 1 2 



3 \2 J V3A Ji-A 
1 8 



(45) 
(46) 



2J2 - A 3( Ji + J2 - A) 



Setting 2te + /ie = wc obtain the critical field A'^^ up 
to order 0{J^,n^ / Ji_) 



A+ - 3 



12Jj_' 



(47) 



where we neglected corrections of order J2 and ft^/Ji. 



4- Phase diagram in the (A, Q,) plane 



We started our analysis in scction lTl A II with the classi- 
cal limit of the Hamiltonian ([7]), J± = ft = 0, where two 
transitions as a function of A were identified. A finite 
hopping Jj_ replaces each of these two classical transi- 
tions with two extended gapless phases (see Fig. [T]) where 
the fluctuations of the magnetization are described by the 
Luttinger liquid theory. As a consequence, at = and 
Jl > one crosses four Lifshitz transitions as a func- 
tion of increasing A. This picture is, however, modified 
if a finite f2 is taken into account as shown in Fig. [21 
which summarizes the phase diagram in the (A, 57) plane 
at some small but finite J±. 

The A axis in Fig. [2] corresponds to a horizontal 
cut through the phase diagram of Fig. [T] The Lut- 
tinger liquid phases (blue shaded) that are generated by 
the finite hopping J± and extend for f2 = between 



A~<A< A+ with i = l,2, are now quenched by the fi- 
nite driving > to only small regions close to A^j^ and 
A^2- (The extension of the Luttinger phases is exagger- 
ated in the figure in order to visualize them) . Both small 
regions surround a commensurate phase with a broken 
Z3 symmetry which is accompanied by a magnetization 
plateau at m = — i. The C-IC transition between this 
commensurate phase and the incommensurate Luttinger 
liquids is of Lifshitz-type, i.e., they are in the standard C- 
IC universality class^^^ so that the magnetization m(A) 
close to these transitions changes with the characteristic 
square root cusps. The values of the Luttinger parame- 
ter at these transitions is K = i; see Eqs. (|22p and (|42p . 
We evaluated the phase boundaries A^^(ri) and A^2(^) 
in lowest order in V, and J±; see Eqs. ((35|) and (|47l). re- 
spectively. The Luttinger liquids themselves are bounded 
by two Kosterlitz-Thouless (KT) transitions into gapped 
incommensurate phases. At these transitions, the Lut- 
tinger parameter reaches the value K = ^ at which 
the operator Sf driving the transitions becomes relevant. 
The positions A^"'"(0) of the two KT transitions, i = 1,2, 
in the limit f2 — >■ were evaluated in Eqs. p5|) and (|43p . 
From the RG flow of the KT transition it follows that the 
phase boundaries start linearly, Af^{n) — A^'^(O) ~ fi. 

Wc did not calculate the phase diagram for larger 
values of fi. From previous works treating similar 
problems,—"— the presence of a special point indicated 
by P3 in Fig.[5]is known where a direct transition between 
the broken Z3 phase and the gapped IC phase takes place 
which is in the Z3 Potts universality class. We assumed 
here that the two Luttinger liquids meet at P3.— In the 
limit Jl ^ J2, the position of this quantum phase tran- 
sition can be estimated to be at 57 ~ A '--^ J2'^ 

Finally, at larger A the gapped IC phase is separated 
from the Z2 broken phase by an Ising transition (I) whose 
boundary A^(ri) was evaluated in Eq. (|34l) . 

We considered the phase diagram for a finite J± and 
evaluated the phase boundaries perturbatively in 51. In 
the limit J± — >• 0, the phase boundaries A~2, A^'^, and 
A+2 wiU merge for ^- 0; see Eqs. ^ and dSJ). It 
was demonstrated explicitly in Ref. that in the limit, 
J± — > 0, the Luttinger liquid phase bounded by A~2 and 
A^"^ nevertheless survives for 57 > 0. Note that the 
evaluation of the phase boundary A^^^ (57) in Eq. (|47p is 
only controlled if 57 <C J_l and the limits 57 — > and 
Jj^ — > do not commute. 



B. Anisotropic spin chain with long-range 
interactions 

In the last section, we demonstrated that the presence 
of an additional next-neighbor interaction gives rise to 
an additional broken Z3 phase that is surrounded by a 
floating Luttinger liquid phase. Now we consider the ad- 
ditional modifications caused by including longer range 
interactions. We can anticipate that an interaction be- 
tween further distant neighbors stabilizes commensurate 
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phases of higher periodicities p breaking Zp symmetries 
of the lattice. 

To prepare the basis for the following discussions, we 
first generalize the treatments of sections fll A 2l and lII A 31 
to classical transition points between two arbitrary ad- 
jacent commensurate phases and study the influence of 
a finite hopping J±. Afterwards we discuss the resulting 
modifications to the phase diagram. 



1. General transitions between commensurate phases at 
finite Jx > and = 

If a Hamiltonian with a finite-range interaction ([6]) 
is considered with n > 2, higher order commensurate 
phases will be stabilized on the classical level. Here, fol- 
lowing closely section III Al we consider the influence of a 
finite hopping Jj_ on the classical transition between two 
arbitrary adjacent commensurate phases with periodicity 
p and p' and magnetization rup and rup' , respectively. 

We assume nip < rripi , and we will also use the relation 
between magnetizations and periodicities 



mp'pp 



1 



nippp 



(48) 



The latter can be derived from results of Bak and 
Bruinsma4i They demonstrated that as a function of 
increasing A commensurate phases are successively sta- 
bilized such that their fraction q/p of up-spins and, 
thus, their magnetization rUp = q/p — 1/2 increases 
monotonously. It is elementary to show that for a given 
commensurate phase with a fraction q/p, the adjacent 
phase with larger fraction q' /p' of up-spins can be found 
by solving the equation qp' = pq' — 1, and choosing the 
solution with largest p' . The choice for p' is only bounded 
by the range of the interaction, i.e., hy p' < n + 1 with n 
of Eq. ([6|) . Substituting the magnetizations for q and q' 
one arrives at the above equation. 

Near the classical transition we can limit ourselves 
to the degenerate states with the successive occur- 
rence of unit cells of the p and p' states only, e.g. 

For a definition of these unit cells 
corresponding to arbitrary fractions q/p and q' /p'-, and 
for a proof that the above limited set of states becomes 
the degenerate ground state manifold at the classical 
transitions; see appendix |B] We again introduce ficti- 
tious particles living on a fictitious lattice and associate 
the unit cell Up with an empty site and the unit cell Upi 
with an occupied site so that the above sequence of unit 
cells is identified with ...001011... The number of par- 
ticles is denoted as N and the number of empty sites 
by Ne, so that L = pN^ + p'N with the length of the 
chain L. The maximum number of fictitious particles is 
< < L/p' and the length of the fictitious lattice is 
Lf = Ne+ N = {L - N{p' - p))/p. The magnetization 
per site is a function of N and given by m = mp+N/ (pL) 
ranging between mp < m < uipi . 

At the classical transition, all states with < A^ < 
L/p' are degenerate. We can evaluate the resulting en- 



tropy S = log fls where is the total number of avail- 
able states. 



L/p' 

E 



L-jV(p'~p) \, 



P 



N-- 



L~N{p'-p) 



N]\ 



(49) 



where each term in the sum is just the number of ar- 
rangements of N particles on Lf = {L — N{p' — p))/p 
sites. 

A finite J± allows the particles to propagate on the 
fictitious lattice. Note that adjacent unit cells with dif- 
ferent periodicity, Up and Upi, can be exchanged by ap- 
plying once the hopping operator with amplitude J_l; see 
appendix [BJ So we obtain the effective Hamiltonian de- 
scribing the propagation of fictitious particles in terms of 
a tight-binding model for hard-core bosons. 



Heff = 



i=l 



J± 

2 



a] a, 



A, 



a] a,; 



(50) 



The chemical potential is obtained with the help of the 
relation between the number of particles A'^ and the mag- 
netization m. The resulting magnetic energy per site 
reads 

e(m) = — m(A — Ac) H — —{m{p — p') + mpip' — nipp) 



7r(m — mpi)p' 



m{p — p') + nipip' — nipp 



(51) 



that has to be minimized with respect to m, de/dm = 
with nip < m < rripi . The finite J± gives rise to a 
compressible phase for fields A in the range 



Ac - pJ± < A < Ac + p'J±. 



(52) 



The magnetization varies continuously within this range, 
and at both boundaries a Lifshitz transition takes place. 
As before wc can evaluate the susceptibility x and with 
the help of the velocity of the domain walls we derive the 
Luttinger parameter within the compressible phase 



K = {{nip — m)p — {nipr — 7n)p') 



(53) 



With the help of the equation Eq. (|3S]), we obtain the 
value of the Luttinger parameter close to the Lifshitz 
transition into a state with periodicity p 



(54) 



This value validates that the commensurate phases are 
stabilized by Umklapp scattering of p spinonsi^ 



2. Higher order commensurate phases for longer range 
interactions 

If a longer range interaction is considered, higher order 
commensurate phases appear on the classical leveli^ For 
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J± 

K = 1 



n 







/ 





















A,l A,2 A,; 



Arf A,5 A 



FIG. 3: Upper panel shows the magnetization as a function 
of A of the classical auxiliary Hamiltonian with a finite range 
interaction J^^^ of Eq. ©. There are five classical transi- 
tions; see Eqs. (|55|l . The spin configurations in the lower 
panel illustrate the structure of the unit cells of the periodic 
commensurate spin states. A finite J± gives rise to Luttinger 
liquid phases (shaded regions) [see Eq. (|52|) ]: the value of the 
parameter K at their edges is given by Eq. (|54p and shown in 
the plot. The figure is not to scale; the extension of phases 
with larger unit cells is exaggerated. 



example, with the interaction J^^^ of Eq. ([6]) with n ~ A, 
there are akeady six classical ground states at J_l = = 
in the vicinity of A = 0. The positions of the phase 
transitions between those states are given hy^ 



Ael = 0, A,2 
Ac3 - 4J3 - 3 J4, A,4 
Ac5 = 3J2 — 2 J3 + 5J4. 



5J4, 

3J2 - 2J3, 



(55) 



The sequence of these classical ground states is illustrated 
in Fig. [31 The figure is not to scale; we exaggerated the 
extension of phases with larger periodicities in order to 
visualize them. For A < Ad = 0, the fully polarized 
state is the ground state with magnetization m = — i. 
For positive A commensurate phases with various peri- 
odicities p arc stabilized with the concomitant breaking 
of a Zp symmetry. There are single phases with peri- 
odicity p = 2, p = 3 and p = 4 with magnetization 
m = 0, — |, — respectively, and two phases with p = h 
and magnetizations m = —j^ and m ~ —j^. The phases 
are arranged such that the magnetization as a function 
of A monotonically increases in a stepwise fashion. Each 
classical transition possesses a macroscopic degeneracy 
given by Eq. (|49l) and thus a finite entropy that will be 
quenched by a finite hopping Jj^. Furthermore, for a fi- 
nite Jj_ as shown in Fig.[3]each classical transition serves 
as the origin of a fan containing a Luttinger liquid phase 
with an extension given by Eq. ([5^ . The edges of each 
fan are Lifshitz transitions where the Luttinger liquid pa- 
rameter assumes the value given in Eq. (j54p . 

It is clear that for sufficiently large J± Luttinger liq- 
uid phases attributed to adjacent classical transitions will 
start to overlap (not shown). This has important conse- 
quence for the full Hamiltonian ^ with the long-range 
interaction J\j.\. As the extension of a commensurate 
state with periodicity p on the = axis is on the or- 
der of 5A - Jp_i + 0( Jp) for Jj_ = 0^ it follows from 
Eq. ((5^ that for Jp_i < Jj^ all higher commensurate 



FIG. 4: Topology of the phase diagram for the Hamiltonian 
Q for a hopping J± > J5 with commensurate phases break- 
ing Z2, Z3, Z4 and Z5 (two smallest lobes) symmetries. Com- 
mensurate phases with higher periodicities are washed out. 
The positions of the critical fields Ad are approximately given 
by Eqs. (|55p . The figure is not to scale; the extension of the 
commensurate phases with larger unit cells is exaggerated to 
make them visible. 



phases of period p will be washed out. This means that 
if, e.g., J5 < J±, one crosses at = as a function of 
A only 10 Lifshitz transitions as in Fig. |3] even if the full 
long-range interaction J^r\ of Eq. ([2|) is considered. This 
simplifies the topology of the phase diagram in the (A, fi) 
plane of the full Hamiltonian ([3|) considerably. 



3. Phase diagram 

The phase diagram of the Hamiltonian ([3]) with the 
long-range interaction J|r| in the (A, fl) plane defined 
by the Rabi frequency Q and detuning A is shown in 
Fig. m For simplicity and illustrative purposes, we as- 
sume a finite hopping of Rydberg excitations on the order 
of J± ^ J5 so that commensurate phases with periodic- 
ity p > 5 are not realized. Note that in typical optical 
lattice experiments the super-exchange interaction is usu- 
ally much smaller J± ^ Jp, with p reaching up to 10'^ for 
deep optical lattices; see discussion in section HVBI The 
restriction to J_l ~ Jp with p = b streamlines the dis- 
cussion and captures all essential features with obvious 
generalizations to larger p. Apart from a broken Z2 Ising 
phase, there arc broken Z3, Z4 and Z5 phases for such 
parameters whose corresponding spin configuration was 
illustrated already in Fig. [31 Up to corrections of order 
J5, the critical fields Ac; are given by Eq. ([55]). 

The phases on the axis £7 = follow from a vertical 
cut through a phase diagram like that of the lower panel 
of Fig. [3| at a finite Jl ^ Jbi so that only 10 Lifshitz 
transitions are crossed giving rise to five Luttinger liquid 
phases (blue shaded) on the = axis. Repeating sim- 
ilar arguments as in section III Al about the stability of 
these Luttinger liquid phases with respect to the Sf op- 
erator associated with a finite f2, we find that the broken 
Z2 phase is bounded by an Ising transition line denoted 
by /. The other broken Zp phases are surrounded by a 
fioating Luttinger liquid phase that extends to finite VL 
and is itself bounded by a Kosterlitz-Thouless transition 
(KT). At the latter transition, the Luttinger liquid pa- 
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rameter assumes the universal value K = 1/8. We will 
discuss the physics and the signatures of this KT transi- 
tion in the next section in more detail. 

The topology of the phase diagram for larger D, follows 
from known results.— i^i^"— First, the KT transition 
line close to Ad connects with the one close to the Ising 
transition so that the Luttinger liquid phase houses all 
gapped commensurate phase of periodicity p > 2. Each 
broken Zp phase possesses a special multicritical tran- 
sition point where the magnetization inside and outside 
the symmetry broken phase is accidentally commensu- 
rate. Whereas for the and Z4 phases this point con- 
nects directly with the polarized phase for large $7, it 
is located within the Luttinger liquid phase for the 
phase (and also for higher commensurate phases that ap- 
pear for smaller Jj_). For periodicity p = 3,4 this mul- 
ticritical point is of Potts type^^i^'^ and for the higher 
commensurate phases it is in the KT universality class 
with K = 2/p^, i.e. the Luttinger parameter at the mul- 
ticritical point is larger by a factor of two compared to 
the corresponding Lifshitz transition; see Eq. ([M)) . The 
transitions from the commensurate phases into the Lut- 
tinger liquid away from these special points are instead 
Lifshitz transitions 4^1^ The width SA and height Sft of 
the dome housing the Luttinger liquid phases can be es- 
timated to be on the order of the next-nearest neighbor 
coupling, 6Ar^ J2 + ©(^3). 



III. DISLOCATION-UNBINDING TRANSITION 
OF THE RYDBERG CRYSTAL 



We now concentrate on the parameter regime of the 
phase diagram where both Rabi frequency, fi, and de- 
tuning, A, arc small compared to the hopping of Ryd- 
berg excitations, J± > fl, |A|. We have seen in the last 
section that for increasing A a first commensurate state 
is stabilized for £7 = only at A Jj^; see Eq. ([5^ . We 
therefore expect that in the regime A <C Jl a contin- 
uum description that neglects any commensurability ef- 
fects should be possible. Such a description is at the focus 
of this section. This allows us, in particular, to elucidate 
the nature of the floating Luttinger liquid phase and, in 
particular, the Kosterlitz-Thouless transition with Lut- 
tinger parameter K = 1/8. 

We start with representing the spin operators with the 



standard Jordan- Wigner fermions, {//, fj} 



Stifle 



J't>, 



— TT-i 2 ~ fi fi 2' 



1 1 



(56) 
(57) 



so that the Hamiltonian ^ becomes 



(58) 



We will again apply the strategy that we analyze first 
the Hamiltonian for $7 = 0, and then consider the effect 
of a finite Rabi frequency perturbatively. For ft ~ 
0, we can distinguish between two regimes: whereas for 
A < A* , with A* to be determined below, the Rydberg 
excitations behave like a dilute weakly interacting gas, 
they form an incommensurate Rydberg crystal for A* < 
A < J±. 



A. Dilute weakly-interacting Rydberg gas 



Consider the Hamiltonian (|58l) for 17 = at very small 
densities of Rydberg excitations (hi) <C 1. In this limit, 
there are only few Rydberg excitations that hardly in- 
teract. Neglecting this residual interaction, the dilute 
gas of Rydberg excitations can be approximated by the 
single-particle Hamiltonian 



dxfix) 



2m 



M fix), 



(59) 



where we have taken the continuum limit of small optical 
lattice constant, a — 0, with fi — > f{x)/y/a, and h = 
P f . The mass m and the chemical potential fi simply 
follow from the single-particle properties of Eq. ((58)) . 



1 

m 



— = Jj_a2, 



^l = A + Jx 



(60) 



As the detuning A increases from negative values, the 
excitation gap, — /i, for a single Rydberg excitation de- 
creases and vanishes at 



A- = -J^. 



(61) 



At A = A^, there is a Lifshitz transition, the Rydberg 
excitations start to condense and form a Fermi sea. The 
critical value (j6ip is in agreement with the results ob- 
tained in the last section; see Eq. (|52p . The energy den- 
sity as a function of the density of Rydberg excitations 
np> is given by the standard free Fermi gas expression 



6m 



(62) 



For /X > 0, minimization with respect to the density gives 
ni^ = \/2m^/TT and, thus, the ground state energy of the 
Rydberg gas 



£rg 



2V2 
"3^ 



m/i 



3/2 



(63) 
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The residual interaction among Rydberg excitations ne- 
glected in Eq. ([5^ can in principle be extracted from 
the two-particle sector of the lattice Hamiltonian ([58| . 
In the low-energy limit, we expect it to have the form 
f\x){dx.f\x))f{x){dxf{x)), where the two gradients are 
required by the Pauli principle. This gives rise to a cor- 
rection to Eq. (p^ on the order of C('n-|j) that is sub- 
leading in the low-density limit nn — > 0. 

For small but finite densities, > 0, i.e., A > Aj, 
the small energy excitations on top of the Fermi sea are 
governed by the Luttinger liquid theory [see Eq. ([20)) ] . 
characterized by the velocity v and the Luttinger pa- 
rameter K. They obey here the relation v = vp/K 
(which follows form the emergent Galilean symmetry), 
where vp = ■\/2/i/m is the single-particle velocity, and 
the Luttinger parameter derives from the compressibility 
d'^e{nii)/dn\ = vk/ K yielding 



(64) 



The Luttinger parameter attains the value K = \ close to 
the Lifshitz transition, — >■ 0, i.e., — > 0"*" in agreement 
with Eq. (j42|) or, more generally, with Eq. (|54| . 

A Luttinger liquid with such a large parameter K be- 
comes immediately unstable in the presence of a S^' per- 
turbation, that is associated with the transverse field fi. 
The operator S*^ is relevant for > 1/8 [see Eq. ([M]) ]. 
and thus opens a gap. So we obtain the result that in 
the small density limit, the gas of Rydberg excitations is 
only gapless for $7 = and gapped for any finite fi. We 
now turn to the situation of larger densities. 



B. Incommensurate Rydberg crystal 

If at = the detuning A is further increased, the 
density of Rydberg excitations increases and the inter- 
action becomes more and more important. If A is suffi- 
ciently large, A > A*, the interaction J\j.\ finally cannot 
be treated perturbatively anymore. In order to determine 
the crossover value A*, we start now from the opposite 
limit and treat perturbatively the kinetic energy of Ryd- 
berg excitations, i.e., the coefficient 1/m in Eq. ([59)1 . 



1. Rydberg crystal for 1/m = Jl = 

Consider the classical limit of the Hamiltonian ([58| . 
In the continuum limit of small optical lattice constant, 
a — >■ 0, it reads 

1 



(65) 

with /i given in Eqs. (|60p . For fi > there is a finite num- 
ber of particles, i.e., Rydberg excitations in the ground- 
state that interact via the long-range interaction J\r\ ■ In 



the thermodynamic limit L — cx), an eigenstate with a 
fixed density nn of Rydberg excitations is 



(66) 



] = ~oo 



where we abbreviated S+{x) = /t(x)e-*'^(^). The state 
I^'rc) depends on the positions, Xj, that we assume to 
be ordered such that ... < Xj-i < Xj < Xj+i < .... 

The state |4'rc) is an eigenstate of the density operator 
n{x) with an eigenvalue given by the Dirac comb 



^(a;)|*Rc> 



OO 

E 

j = -(X> 



S{x-x,)\^nc)- (67) 



2. Ground state energy of the Rydberg crystal 



In the ground state I^'rc) the Rydberg excitations are 
equally spaced so that the positions Xj of the wavcfunc- 



tion (1661) assume the values x 



(0) 



I* 



(0)n 
RCi 



- n s^i'^rm 



(68) 



The corresponding energy density e(nfl)|4'^Q) = 
X^ciassj^Rc) is a function of nn and reads 

^i^R) = I > , J\,n\/(nna) - A* | "fl- 




(69) 



Using Eq. ^ for the interaction, the energy simplifies to 

e(nfl) = jR({a){nRa)°'nR - ^iur, (70) 

that has to be minimized in order to determine ur, 
de/driR = 0, 

l/a 



1 



"^=aU.(« + l)C(a)j ' ^''^ 

with C the Ricmann zeta function. Taking e{nR) at its 
minimum one finally obtains the ground-state energy for 
the Rydberg crystal state, 

a 

£RC = -- — TTTTTTr — Mi/^ -1/^ . (72) 



(a + l)[(l + a)C(a)]V« jV«« ' 

that is non-pcrturbative in the interaction Jr. 

Comparing the energies of the Rydberg gas, Eq. ([63]) . 
and the Rydberg crystal, we obtain for the crossover 
value of the detuning 



A* 



A^ , with /i* ~ Jl 



Jr 



(73) 



and A^ = -J±. Note that A* < for J± < Jr and, 
consequently, the incommensurate Rydberg crystal exists 
in a finite regime A* < A < Jj_ between the Rydberg 
gas at smaller and the commensurate Rydberg crystal at 
larger detuning A. 
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3. Degenerate perturbation theory in the kinetic term 

The incommensurate Rydberg-crystal ground state of 
Eq. breaks the translational symmetry of the contin- 
uum Hamiltonian. Consequently, quantum fluctuations 
induced by a finite hopping J± will give rise to phonon 
excitations. 

Besides the simple shift in the expression for the chem- 
ical potential, fi, a finite J± results in a kinetic energy 
with mass m as defined in Eq. (|60| . Consider the crys- 
talline wavefunction in Eq. (|66|) where the positions of 
the excitations Xj differ from their equilibrium positions 

by Uj ~ Xj — x^^"^ . The eigenvalue problem in the pres- 
ence of the kinetic term, ["Hciass + / dxp {—dl/{2'm))f — 
^'ll^Rc) ~ 0, then reduces to finding the solution of the 
first- quantized Hamiltonian 



(74) 



where the momentum and position obey the canonical 
commutation relation [xj,pi] = iSij, and Ni^^ = nuL. 

Expanding the Hamiltonian in second order in the 
small deviations Uj we obtain H w Le{nj^) + -ffphonon 
with the energy £{nii) of Eq. ([70| and the phonon Hamil- 
tonian 



phonon 



j 



1 



(75) 



where J|"| = d'^ J^r\ / d\r\'^ . In order to obtain the asymp- 
totics at lowest energies, we consider a second contin- 
uum limit, i.e., we confine ourselves to length scales much 
larger than the mean distance between Rydberg excita- 
tions, 1/jT.fl. In this limit, we can use Ui — >■ u{x) and 
Pi ~^ p{x)j with the commutation relation [u{x),p{x')] ~ 
iS{x — x')/nji, and the phonon Hamiltonian is approxi- 
mated by 



phonon ~ '^R 



dx 



1 B 



2m 2 (anfi)' 



(76) 



where we kept only the lowest derivatives of the u(x) field 
and we abbreviated 



^ ■J\'m\/{nRa''"T- 



2 '^i'^ + ^)JR'm^ 



m — 1 



i\m\/{nna)r+^ 
a(a + 1) Jfl(n,ja)2+"C(a). 



(77) 



Introducing the dimensionless fields u{x) = 
~(j){x)/ (irnj^) and p{x) = —dj;0{x) with the canon- 
ical commutation relations [-(f>{x),dx'd{x')] = iS{x — x'), 
the phonon Hamiltonian |76p can be recast in the 
standard Luttinger liquid form of Eq. (|20|) with velocity 
V and Luttinger parameter K obeying 

Kv=^, ,2 = 1^. (78) 





FIG. 5: Worldlines of the incommensurate Rydberg crystal. 
Vertical axis shows position x in units of 1/ ur and horizontal 
axis shows time t. The wiggles indicate phonon excitations. 
The left panel shows a single dislocation and the right panel 
a dislocation pair. The unbinding of such pairs leads to a 
melting of the crystal at the KT transition. 



Using the explicit expressions for B, for the mass m given 
in Eq. (pO| . and Eq. ((7T|) for Ufj we obtain 




(79) 
(80) 



The Luttinger parameter is small. A' ^ 1, in the incom- 
mensurate Rydberg crystal regime /i* <C fJ-. It increases 
and becomes of order one if the chemical potential is low- 
ered down to the crossover value /i* of Eq. ([73|) . and it 
thus smoothly connects with Eq. (|64p of the dilute Ryd- 
berg gas regime. 

As a consistency check, we might consider the correc- 
tion to the ground state energy due to the phonon excita- 
tions. This is of order Se ^ v J dkk where the momentum 
integral is cutoff at n b, so that one obtains for the relative 
correction 




(81) 



Note that the relative correction is on the order of the 
Luttinger parameter Se/eRc ~ K. So the correction is of 
order one and the perturbative treatment in the collective 
variables u,; and pi breaks down if the chemical potential 
is lowered down to the crossover value /i* of Eq. ([73]) . 



^. Stability analysis with respect to f2 perturbations: 
dislocation-mediated melting transition 

In the presence of a finite Rabi frequency Ryd- 
berg excitations can be added to or removed from the 
crystalline arrangement. In the limit a — > 0, with 
Ji — > /(x')/y^, this perturbation can be written as 



5% 



dx{S+{x) + S-{x)) 



(82) 



where S+{x) = /t(x)e~*'^(^) and S~ {x) = {S+{x))^ . In 
the thermodynamic limit, the operation of S^{x) on the 
Rydberg-crystal wavefunction, Eq. (j66p . i.e., adding a 
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Rydberg excitation can be expressed in terms of a first- 
quantized operator (neglecting Klein factors) , 

S+{x) ~ : e'(^-^")P"+5:?=ioo*(^.+i-=^.)p. (33) 

where the normal ordering, : O :, is defined such that 
the momentum operators should be placed to the right 
of all position operators. The index n for a given x is 
defined such that x„ < x < x„+i. As we require that the 
particle positions are ordered, the addition of an excita- 
tion at position x involves a collective shift of all particles 
with positions Xj < x, that is just realized by the sum 
of momentum operators in the exponent. The above ex- 
pression simplifies in the continuum limit, 

S+{x) ~ _^e*^-- = ^^e-'^("), (84) 

where we used the continuum version of the momentum 
operator and its relation to the B field as defined in the 
last paragraph. We have also set 6'(— 00) = 0. The per- 
turbation ([5^ on the Rydberg crystal thus corresponds 
to the first-quantized operator 

5H = -^=^ [ dxcos{e), (85) 
V27ra J 

[see also Eq. ([M]) ]. This operator can create or destroy 
dislocations in the incommensurate Rydberg crystal. 

Fig. [5] shows the worldlines of Rydberg excitations. In 
the ground state these worldlines are on average equally 
spaced with a distance given by l/rifl. The wiggles 
around the equilibrium positions illustrate the small- 
amplitude, phonon excitations. If the operator S~ is 
applied on the Rydberg crystal, a Rydberg excitation is 
destroyed at a certain time and a dislocation is created; 
see left panel of Fig. [5] The remaining excitations will ac- 
commodate themselves at a later time but the averaged 
distance has slightly increased and, as a consequence, the 
resulting state is energetically unfavorable. In order to 
recover the ground-state configuration, a particle has to 
be added at a later time by applying the operator 
giving rise to an anti-dislocation as shown in the right 
panel. If the incommensurate Rydberg crystal is stable, 
the dislocation and anti-dislocation in the worldlines are 
bound and, in contrast to the phonons, correspond to 
large-amplitude excitations. This is the case if the op- 
erator of Eq. (|85p is irrelevant in the RG sense which 
implies K < l/8i^ As the Luttinger parameter K 1 
for fi fi* , the incommensurate Rydberg crystal is sta- 
ble to such small CI perturbations. The dislocation pairs, 
however, unbound for K > 1/8 and, as a consequence, 
the Rydberg crystal melts in a Kostcrlitz-Thoulcss tran- 
sition. We can estimate the position of the KT transition 
using Eq. ([50]) . The Luttinger parameter becomes of or- 
der one for ^ ~ and the location of the transition 
thus coincides with the crossover between the Rydberg 
gas and the incommensurate Rydberg crystal. As we 
argued previously, from the separatrix of the RG flow 



for the KT transitions^ follows that the phase boundary 
starts linearly in i.e., AKT(f^) — Akt(O) ^ ^■ 

IV. SUMMARY AND DISCUSSION 

Combining the results of section HIl and Hill we summa- 
rize in the following the topology of the phase diagram 
for one-dimensional Rydberg atoms, that are governed by 
the Hamiltonian (jS]). We first outline the global structure 
of the phase diagram before we focus on the experimen- 
tally relevant regime and discuss possible experimental 
consequences. 



A. Phase diagram 

The topology of the phase diagram depends qualita- 
tively on the relation of A and il with respect to the 
hopping J_L and, especially, the next-nearest neighbor 
interaction J2 = Jii/2"; see Fig. [SI 

On scales Jj., J2 < |A|,fJ < Jr [see Fig. El^a)], the 
effect of hopping J± and higher order interactions J|r| 
with r > 1 is negligible and the Hamiltonian ([3]) reduces 
effectively to an Ising model in a longitudinal and trans- 
verse magnetic field. Consequently, one obtains a single 
Ising transition located^ at oc A between a paramag- 
net and an antiferromagnet with a Rydberg excitation at 
every second lattice site, which breaks a Z2 symmetry of 
the lattice. For smaller < J2 the Ising transition line, 
however, deviates from the straight line and approaches 
a value on the 17 = axis on the order A ^ J2 [neglecting 
corrections of order ©(Js)]. 

At intermediate scales J± ^ |A|,r2 ^ J2, an addi- 
tional dome close to the Ising transition appears that 
contains higher commensurate phases. Fig. |6ljb) shows 
clearly the presence of a lobe with a broken Z3 phase 
where a Rydberg excitation is present at every third site 
of the optical lattice. The width and height of the lobes 
decreases fast with the order of commensurability so that 
higher-order commensurate phases are barely visible on 
the scale of Fig.lSJb). In fact, most of the commensurate 
phases of high order are not realized due to the finite 
hopping J±. As a rule of thumb, one can say that only 
those commensurate phases breaking a Zp symmetry of 
the optical lattice appear in the (A, H) plane which full- 
fill J]_ < Jp-i; all phases with higher commensurabil- 
ities p' with Jp>-i <C J_l ai'c washed out; see discus- 
sion in section IIIB II In addition, the extension of the 
broken Zp phases on the 17 = axis is approximately 
on the order of 6A ^ Jp-i + 0{Jp), and thus strongly 
decreases with increasing p. The regions between these 
commensurate phases is filled with an incommensurate 
fioating phase. Here, the Rydberg crystal has a period- 
icity that is incommensurate with the underlying optical 
lattice and, in addition, supports low-energy phonon ex- 
citations described by the Luttinger liquid theory. The 
existence of this floating phase was already pointed out 
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FIG. 6: Sketch of the phase diagram in the (A, Q) plane of the 
Hamiltonian ((3| close to the origin but on various scales: (a) 
for |A|,J7> Ja, (b)for Jx < | A|, J7 < Ja and (c) |A|,n < Jx; 
see text. There are gapped commensurate (yellow), gapless 
incommensurate (blue) phases, and a gapped polarized phase 
(white). The commensurate phases break a Zp symmetry of 
the optical lattice and their extension decreases fast with the 
order of commensurability p. The commensurability of the 
last commensurate phase at small A in (c) is determined by 
Jp ~ J±, i.e. p ~ [Jr/ Ji^Y^"' ■ The abbreviations / and KT 
indicate an Ising transition and a Kosterlitz-Thouless tran- 
sition line, respectively. Points P„ denote critical points of 
n-state Potts universality. 



range Ising model is restored. 

The transition between the incommensurate to a cer- 
tain commensurate Rydberg crystal state is in the usual 
universality class of commensurate-incommensurate (C- 
IC) transitions^!^ except at special points of acciden- 
tal commensur abilities. Furthermore, we find that the 
boundary of the dome is given by a Kosterlitz-Thouless 
(KT) transition where the incommensurate Rydberg 
crystal melts due to the proliferation of dislocationsi^i^ 
The Luttinger parameter at this KT transition takes the 
universal value if = 1/8 with the concomitant universal 
correlation tails for Rydberg excitations. There are two 
special multicritical point on the boundary of the dome 
where the magnetization is accidentally commensurate 
with the and Z4 phases denoted by P3 and P4, re- 
spectively. At these two points, the transition is not of 
KT type but instead in the three- and four-state Potts 
universality class42i^ 

Finally, at very small scales |A|, ^ Jj^ [see Fig. [Sic)], 
commensurability effects are negligible and an effective 
continuum description is possible; see section IIIII This 
allows to describe the physics at the lower left edge of 
the dome where the KT transition hits the = axis. 
The floating incommensurate phase can be here captured 
within a standard strong-coupling perturbation theory 
familiar, for example, from the treatment of Wigner crys- 
tals. The phonons of the incommensurate Rydberg crys- 
tal give rise to gapless excitations and the Luttinger liq- 
uid parameter can be determined. We find that the cre- 
ation of dislocations by the Rabi frequency is irrele- 
vant and the floating phase is stable as long as the de- 
tuning exceeds A > A*, where A* = —J±_ + fi* with 
M* ^ J±{J±/ Jr)'^^^""'^'' ■ For smaller detuning, the Lut- 
tinger parameter reaches K = 1/8 and a dislocation- 
mediated (KT) melting of the incommensurate Rydberg 
crystal sets in. If A increases to values of order pj± [see 
Eq. ([5^ ]. where p ^ (Jr/ Jj^y^" follows from the cri- 
terion Jp ^ Jl, the floating phase becomes unstable as 
well and a flrst commensurate state with periodicity p 
takes over. 



in Ref. [2^ This floating phase is microscopically stabi- 
lized by the hopping of Rydberg excitations. Although 
the bare value of the hopping is very small, J± <^ |A|, fJ, 
it is strongly enhanced by renormalization effects gener- 
ated by the laser drive so that the floating phase 
even extends to large values of < J2. This in fact im- 
plies that the phase diagrams in Figs. IHl^a) and (b) even 
apply to the model Eq. ^ in the limit Jl ^ 0, 

nj^^o = E (S- + 1/2) (5| + 1/2) (86) 

j 3 

In this limit the full fractal structure of the classical long 



B. Experimental parameters and signatures of the 
two melting transitions 

In order to estimate the location of a typical ex- 
periment on Rydberg atoms in the phase diagram we 
now discuss typical values for parameters. In case of a 
van dcr Waals interaction between Rydberg atoms the 
largest energy scale Jr is given by Jr = C^/a^ with 
Cg ~ 500GHz /XTO^ for heavy alkali atoms with principal 
quantum number of order n « 70i^ If the interaction 
is instead of dipolar type one has Jr = Vdd/o? where 
Vdd ~ IGHz/xm*^ for heavy alkalis with n « 15i^ Us- 
ing an optical lattice constant a « 0.5/xm one obtains 
Jr « 32THz and Jr sa 8GHz for the van der Waals and 
the dipolar interaction, respectively. Rabi frequencies of 
order Vl « lOOkHz have been achieved^i that is on the or- 
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der of the r^^ nearest neighbor interaction, = Jj^/r"', 
with r « 26 and r « 43 for the two cases, respectively. 
On the other hand, for a deep optical lattice with po- 
tential depth of 20Er we can estimate the hopping to be 
J± « 10^^ Eri^ For this rough estimate we've assumed 
that the on-site interaction energy between Rydberg and 
ground-state atoms is of the same order as the interaction 
energy between two ground-state atoms. For a typical re- 
coil energy Er ~ lOkHz this yields J± w 10~^Hz that is 
very small compared to all other energy scales, justifying 
the approximation of a frozen Rydberg gas. 

Taken together, this would place a typical experiment 
in the parameter regime pictured in Fig. Elb) . An impor- 
tant requirement for the observation of different phases 
in this phase diagram is, however, that the corresponding 
ground state energy is much larger than the decay rate of 
the Rydberg excitations, 1/r. In particular, the life-time 
T should be much larger than the time of several Rabi 
cycles in order to enable the crystallization of Rydberg 
excitations. Even then, it will be very difficult to iden- 
tify experimentally the universality classes of transitions 
between different phases of Fig. |6] because due to critical 
slowing down very long life-times of Rydberg excitations 
would be required, especially, to distinguish between the 
Kosterlitz-Thouless and a possible first order transition. 
Moreover, it is important to note that experiments with 
quasi one-dimensional atomic gases usually operate in a 
configuration where a large number of independent, par- 
allel ID tubes are arranged in a 2D array with the number 
of atoms per tube typically around ~ 100. In order to 
observe a crystalline arrangement of Rydberg excitation, 
one needs to excite a sizable fraction of the atoms, say, 
about 1/10, so that one has to operate in a regime where 
the parameters ^ A ^ Jn/lO" are accessible. For ex- 
perimental protocols for cooling the system in order to 
realize Rydberg crystallization we refer the reader to the 
Iiteraturei22i2^ 

The important finding of this work is that as function 
of increasing detuning A at fixed one approaches first 
a Kosterlitz-Thouless transition to enter the regime of an 
incommensurate Rydberg crystal. In principle, this KT 
melting transition can be detected by Bragg scattering of 
light, which can be achieved by a resonant de-excitation 
of the Rydberg atoms. If the Bragg signal is measured 
along a ID tube axis, it is just proportional to the ID 
static structure factor 



S{q) = / dxe"'^S'{x)S'{0)), 



(87) 



incommensurate Rydberg crystal phase, one obtains 
1 1 



Siq) 



|g + 27rnH|i-2^ \q - 2Trnn\^~2K 



(88) 



where the Luttinger parameter assumes its maximal 
value K = 1/8 at the KT transition. The structure fac- 
tor is thus expected to show a quasi-Bragg peak if the 
wavelength of the light matches the lattice spacing of the 
floating Rydberg crystal. As the quantum fluctuations 
are even further suppressed so that one reaches the C-IC 
transition and a commensurate Rydberg crystals forms 
(yellow shaded in Fig. [6|), true long-range crystalline 
order finally gives rise to well-defined (delta-function) 
Bragg peaks. In this way, the two melting transitions 
could be experimentally detected. 
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Appendix A: Driven Rydberg atoms in the 
superfluid regime 

So far we have discussed the ground state phase dia- 
gram of driven Rydberg atoms in a one-dimensional op- 
tical lattice deep in the Mott insulating regime, where 
the Hamiltonian ([1]) can be mapped to the effective spin 
model Virtual hopping processes, to-, via high-energy 
states involving doubly occupied sites contribute to the 
nearest-neighbor spin-exchange. Evaluation of this con- 
tribution in second order in the hopping and in zeroth 
order in one obtains 



4ttt. 



Jz 



2(t2 



4t2 



n 



n 



u 



(Al) 
(A2) 



where A = ^.Ti/q is the wavelength of the light and the 
spin, = nR — \, is related to the density of Rydberg 
excitations n/f. In the Rydberg gas phase, the struc- 
ture factor does not exhibit any Bragg peaks at some 
finite momentum q. However, as the KT transition is 
reached and one enters the floating crystal phase (blue 
shaded in Fig. [5]), quasi-Bragg peaks will appear. From 
the algebraic decay of correlations in the floating 



In Eq. ((3)) , we neglected the longitudinal part Jz because 
it only gives rise to a negligible correction of the Rydberg- 
Rydberg interaction on nearest neighbor sites, which is 
several orders of magnitude larger. At this order one 
obtains the standard expression for the generated spin 
exchangci^i^^ As the initial and the virtual states have 
the same number of up and down spin, the energy dif- 
ferences appearing in the denominators do not depend 
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on the laser detuning A. The correction to these expres- 
sions at finite originates from exciting and de-exciting 
a Rydberg atom or vice versa while in the virtual state 
it is located on a doubly occupied site. This yields a cor- 
rection on the order of O ^fjl^-j • For > |A|, this cor- 
rection becomes large and will dominate over Eq. lAil). 
However, it turns out that in this parameter regime, the 
spin-exchange J± is strongly renormalized by processes 
within the low-energy Hilbert space so that its bare value 
can in fact be neglected. If two neighboring atoms are 
both excited or de-excited with an effective J± on the 
order of is generated^ that does not involve dou- 

bly occupied sites and, therefore, is much larger than the 
bare one. For our purposes, we can therefore restrict 
ourselves to the lowest order expression (jAip . 

Now we briefly consider the opposite limit, where the 
ground-state atoms arc in the superfluid regime and 
form a (quasi-)condensate which has been addressed in 
a recent experiment on Rydberg atoms in a ID lat- 
tice.— In this case, we can apply the replacement bj^ — > 
y/no + ajie~'^^'^ with the mean density of ground-state 
atoms, no ^ nji, assumed to be much larger than the 
number of excited Rydberg atoms rt^,, and aj^ and 9ji 
account for fluctuations of the density and phase, re- 
spectively. In the ID case, these fluctuations destroy 
true long-range order leaving only a quasi-condensate of 
ground-state atoms with algebraic correlations. The part 
of the Hamiltonian describing the the excited Rydberg 
atoms reads 

^ = - II (^]t^J+it + h-c-) + -^b-ei ^V^n V^T 
j 

-lj2(^-^^ni^a + ^u))^jt^ (A3) 

j 

In writing the above Hamiltonian, we assumed a strong 
on-site repulsion U^^ between excited Rydberg atoms 
which implies the hard-core constraint Uj-^ < 1. If 
we were to neglect the fluctuations of the (quasi)- 
condensate, Oj^ = <Tji — 0, the above Hamiltonian could 
be mapped to the same pseudospin-Hamiltonian of Eq. 
© with J±_ -> and 17 ^no ^ ^^d A A-2t/f;rio- 
Within this approximation, we would expect the same 
ground-state phase diagram as the driven Mott-insulator. 
The main effect of the condensate is a scaling of the Rabi- 
frequency Q, with y^no, as observed experimentally.— In 
general, however, the fluctuations will modify the phase 
diagram. For = 0, we expect the density fluctuations 
CTi to influence the phase diagram only quantitatively and 
not qualitatively at least for incommensurate fillings tiq 
because cr^ couples to an oscillating Rydberg density pro- 
file Uj^. At finite VL, on the other hand, the presence of 
phase fluctuations Oj^ reduce the relevancy of f2 operator 
in the rcnormalization group sense, so that the incom- 
mensurate Rydberg crystal will be stabilized. In addi- 



tion, the Luttinger parameter at the Kosterlitz-Thouless 
melting transition, see Fig. [HI will differ from K = 1/8. 
The influence of fluctuations on the Ising transition of 
Fig. Oi is probably more complicated and left for a fu- 
ture study. 

Appendix B: General particle description for general 
transitions between commensurate crystal phases 

In this appendix we provide a general definition of 
the fictitious particles introduced in Sec lII B ll and show 
that they form a subspace that becomes the degenerate 
ground state manifold at the classical transitions between 
commensurate states. We consider finite range interac- 
tion specified by Eq. ([5]). Before diving into technicalities, 
we outline our main steps and results. From the results 
of Hubbard^ and Bak and Bruinsma^ it follows that 
all states with fraction q/p of up spins with p < n+ \ [n 
is the interaction range in Eq. ^] have a finite stability 
region and occur in increasing order upon increasing the 
longitudinal field A. Given the fraction q/p, in order to 
minimize the interaction Eq. ([6]), the up spins arrange in 
a particular crystalline pattern with periodicity p which 
maximizes their distance. The general formula for the 
unit cell defining this pattern is given^ in Eqs. (|B1|) and 
(|B2p . For example, for the fraction q/p = 2/5 this gives 
Ub =itiit- Next we consider a transition between two 
crystalline phases q/p — >■ q' /p' with p,p' < n -f 1, and 
with unit cells Up and Up'. Consider the set of states 
formed by arbitrary sequences of the two unit cells, e.g., 
...UpUpUpiUpUp.... We will show that this set of states 
is the degenerate ground state manifold at the classical 
q/p ^ q' /p' transition. Furthermore, associating the two 
types of unit cells Up and Upi with either particles or holes 
living on a fictitious lattice (e.g., the state 
denoted as 00100), we calculate the particle hopping am- 
plitude (e.g., the amplitude for 00100 -> 00010) at finite 
J±. The result, leading to Eq. ([50| . is that the hopping 
of the fictitious particles is obtained by means of a single 
operation of J± on the original spins. 

Definition of the unit cell Up in a crystalline 
configuration: Consider the crystalline state with frac- 
tion q/p of up spins. Its unit cell of length p, denoted Up, 
is found using the following prescription)^ (i) define the 
integers fc, no, ni, by 

p/q = no + ro, 
|l/ro| = ni+ri, |l/ri| = n2 + r2, ... 
\l/rk-2\ = nfe_i+rfc-i, \l/rk-i\ = rik, (Bl) 

where for all s, — i < < i. (ii) Define the sequences 
Xi,X2,...,Xk+i and Yi,r2,...,i^fc by 

Xi = no, yi = no + ao, 

= [X,]"'-iy„ y+1 = [X,]"-+"'-iy, (B2) 

where ofj = ^i/kj| ~ il- (iii) Then, the required unit 
ccU is Up = Xk+i- 
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TABLE I: Crystalline states stabilized for a finite range in- 
teraction Eq. (|6]| with n = 13. 



In a sense k measures the degree of "fractallity" of 
the crystalline state: for fc = unit cells have the form 
Xi = riQ, consisting of no — 1 down spins followed by one 
up spin, e.g., Xi = 3 =4.it- For k = I, unit cells have 
the form X2 = [Xi]"i~^Yi, consisting of ni — 1 repeti- 
tions of the string Xi followed by the string Yi = XizLl, 
e.g., [3]3[4] - [3] [3] [3] [4] =;;n;mm;t. similarly, for 
a given k the unit cells X^+i consist of a number of rep- 
etitions of the string Xk followed by the string Yfc. To 
illustrate this construction we display in table U all frac- 
tions and corresponding unit cells for n = 13. 

Degenerate unit cells Up and Upi at classical transitions: 
It is elementary to show that at transitions [between clos- 
est fractions < q/p < q' /p' < 1 with p,p' < n + 1] k 
either (i) remains constant, k = k' , or (ii) k ~ k' = ±1; 



see table HI For transitions of type (i) we can write 

{Up,Up,)^{Xkn,Ykn), (B3) 

where {X^n , Y^n } are given from Eq. (jB2p evaluated for 
a third fraction, q" /p" = {q + q') /{p + p') with degree 
of "fractallity" k" (that has the unit cell Xk"Yh>> which 
is not stabilized for the considered interaction range); 
For transitions of type (ii) we may write {up,Up'} = 
{Xk+i,Xk}, where {Xk+i,Xk} are given from Eq. ([62]) 
evaluated for the (smaller) fraction q/p. [in case (ii), it 
is sufficient to consider the case that A:' = fc — 1] . 

Degeneracy of the sub-Hilbert space: We associate Up 
{upi) with particles (empty sites) on a fictitious lattice 
and set the magnetic field A to the transition between 
the corresponding commensurate crystalline phases. We 
need to show that at Jl = 57 = 0, states with arbitrary 
sequences of Up and Upi , corresponding to arbitrary par- 
ticle states in the fictitious lattice, are degenerate. Since 
the condition of convexity is assumed, upon increasing 
A all sates with rational density of up spins become a 
ground state (possibly degenerate) at some field. If the 
interaction is of infinite range than every fraction ac- 
quires a plateau4i For finite range interaction, however, 
the two crystalline states .. (maximal particle 

occupation) and (zero particle occupa- 

tion) become degenerate at their transition field also with 
an infinite set of crystalline phases, e.g.. 
The unit cells of such periodic arrangement correspond 
via Eqs. (jBip and (|B2I) to fractions in the open interval 
(q/p, q' /p') which might have arbitrary large k values; 
these fractions necessarily have a denominator > 71 + 1 
and hence are not stabilized over a finite range of longi- 
tudinal magnetic field A. Hence, for any density of par- 
ticles there exist degenerate particle states. For example, 
if fc = A;', then we can construct states with arbitrary low 
particle density of the form Up(Mp')" with any integer n, 
or states with arbitrarily large particle density — >■ 1 of the 
form {up)"up'. However, those states are very special as 
they correspond to some periodic crystalline arrangement 
of the spins. 

To show that all particle configurations are degener- 
ate, it remains to be shown that one can start with 
one of those periodic states and move particles around 
without any energy change. The basic such process 
that we need to consider is the single particle hopping, 
. , .Upf Up , . . y . . .UpUp' .... We will consider transitions of 
type (i) between two states with the same value of k 
where Up = Xk" and Up' = Yfc" with k" = k + 1 = k' + 1; 
see Eq. (jB3p . Case (ii) can be addressed along the same 
lines. We start with the initial configuration 

|z) - ...rfc.+iXfc+i. .. = ... 01..., (B4) 

which contains a Yfc+i string followed by an X^+i string, 
surrounded by a particular configuration of either one of 
those strings, hidden in the (...) which is fixed to define 
the state, but can be arbitrary, and we would like to 
arrive at the final state 

|/) = ...Xfc+iyfc+i... = ...10..., (B5) 
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which coincides with the initial state except for the ex- 
change Yfc+iXfc+i — >■ Xfc+iYfc+i. Our purpose below, is 
to show that this exchange process consists of a single 
operation of and that it does not cost energy, inde- 
pendent of the unspecified (...) configurations. 

We use Eq. (|B2|) to express Xk+i and Yfe+i in terms of 
Xk or Yfe. We obtain 



\i) = ...[Xkr^YkXk[Xkr>'Yk..., 
\f)^...[x,r''XkYk[Xkr''Y„..., 

for ak = —1, or 

\i) = ...[Xkr''XkYk[Xkr>'Yk..., 
If) ^ ...[Xkr-YkXk[Xkr''Y,..., 

for ak = 1, where = 71^ — 1-1- 



(B6) 



(B7) 



= min{7ii — 

— 1 -I- ai}. Using Eq. (jB2p we continue successively 
decomposing Xk and Yfc into and Yfc-i until we 

reach the Xi and Fi unit cells. The result is 



\i) = ...BY^XiA..., 
I/) = ...BX.YiA..., 



for nLi(- 



= 1, or 

\^)=. 



..BXiYiA..., 
..BYiXiA..., 



(B8) 



(B9) 



for niLi(~"j) = 1' where 

k 

A = \{[x,rm, 

i=l 
k-1 



B 



(BIO) 



Furthermore, YiXi = [;]"'"(tl)[l]™" t if = -1, or 
Y^Xi = [4.]™''(it)[i]"'" t if «o = 1, hence we obtain 



for nto(- 



|j) = ...i3[;]""(n)[;]™n A.., 
I/) = ...i?[;]""-'(;t)[;]'"n^-, 

ttj) = 1, or 

|j)-...B[i]™"(it)[i]'"nA.., 
I/) ==...i3[i]™«(n)[i]"nA.., 



(Bll) 



(B12) 



for ni=o(-"») = 1- From Eqs. ((BTTI) and ((BT2]) we see 
that 



{^\n\f) = {f\m^-^-^, 



(B13) 



where is given in Eq. ([3]), namely, the hopping ampli- 
tude is just that of flipping the two spins in parentheses, 
achieved by a single operation of J_l. We proceed with 
the proof that the interaction energy of the initial and 
final states are equal. 



The unspecified (...) to the right in Eqs. (|B4[) and 
(|B5)) . either starts with Xk+i = [XkT>'-'^Yk or Yu+i = 
[Xfe]"^-^^+"^-yfc, therefore it begins with [Xfc]™^ Simi- 
larly, the undetermined (...) to the left must start with 
[XkY^'^Yk- Thus, we can display more spins for the same 
state given in Eq. (|B6p . 



ii) = ...[Xkr^Yk[Xkr'{XuYu) 
[Xkr^Yk[XkrK.. 



(B14) 



(for ak ~ ^1)- We may therefore successively intro- 
duce in the right side (...) in Eqs. (jBlip and (|B12|) the 



combinaiton [Xk-i]"^''^^ , then [Xfc_2]' 



till [Xi 



Similarly, we successively add in the (...) to the left 
the strings [Xfc_i]"'^'-in_i, then [Xfe_2]"'^-^Ffe-2, till 
[Xi]™iYi. The result, for ULii^^i) = 1' is 



|z) = ...^i?[i]""(n)[i]™n^S..., 

\f) = ...AB[iroamr°tAB.... 



(B15) 



The unspecified (...) to the right in Eqs. (jBlsp . either 
starts with Xi = [i]"""^ f or Fi = [|]»o-i+ao ^ence 
it begins with = H ... i- Simflarly, the unspecified 

mo 

(...) to the left must end with t- Hence we obtain 



I*) = ...c{n)c..., 

I/) = ...C{it)C..., 



(B16) 



for n 



1=0 V 



1, or 



\i) = ...c(;t)c.... 
If) = ...c(n)c..., 



for n 



i=0\ 



-ai) = 1. Here, 

c = [;]'"" t AB[iY 



(B17) 



(B18) 



What we have done in Eqs. ()B16P and (|B17P is to maxi- 
mally determine the spins to the left and right of the ex- 
changed up and down spins shown in parentheses, which 
physically follows from the crystalline underlying inter- 
action. The interaction of the up spin in parentheses 
with all other up spins within the string displayed in 
Eqs. (jB16p and (jB17[) is equal for the final and initial 
states. This follows because the string C is related to it- 
self by left-right inversion. This property, C = C, where 
X is the left-right inversion of a string X, follows from 
the induction proven set of identities 



[t]Y[Y,x, = x,Y[%^M 



(B19) 



(B20) 
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for any i > 1. Furthermore, due to the finite interac- 
tion range, J^^^ = for \r\ > n, the up spin in paren- 
theses does not interact with any spin out of the dis- 
played sequence in Eqs. (|Blip and (|B12p . as we now 
demonstrate. Denoting the total number of spins (up 
or down) in a string X by Lx, we have \r\ > Lc + 1- 
Furthermore, Lc = L^Xk+i]lYk+i] ~ 2, which follows by 



comparing Eqs. (jBlip . (|B12p . and (|B18|) . Since a crys- 
talline state with unit cell [^fc+i] [^fc+i], corresponding 
to fraction {q + q')/{p + p'), is not stabilized, we have 
L[Xk+i]lYk+i] > n + 2. This implies |r| > n + 1 hence 

(n) 

J|^l = 0, completing our proof that the interaction en- 
ergy in the initial and final states are equal. 
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